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I. INTRODUCTION 

Despite the fact that the QCD Lagrangian can be 
written down in a line, deriving the huge plethora 
of strongly interacting phenomena that are seen in 
accelerators (and in everyday life) has remained a major 
challenge for over half a century. Especially at low 
energies, QCD still remains intractable by traditional 
perturbative methods. A new window into QCD physics 
has recently opened after the experiments carried out at 
RHIC during the last decade.'^ Now, the new collision 
experiments at RHIC and LHC (see for example Ref. [2] 
for early results) will allow obtaining much more precise 
informations about the phase diagram of QCD in the 
region of high temperature and small baryon density. 
During the last years, evidence of the formation of a 
quark-gluon plasma has been growing, in the form of a 



2 



transient object called the "fireball" , which expands and 
eventually explodes into hadrons.-^- 

String theory, which was originally devised to explain 
the low energy phenomenology of QCD, has ended 
up producing one of the most valuable theoretical 
methods to explore strongly coupled systems: the 
AdS /CFT duality. While its prototypical example, 
M — 4: supersymmetric Yang-Mills theory, is far from 
being a realistic analog of QCD, its finite temperature 
counterpart shares with QCD many properties such as 
Debye screening, a deconfined phase, etc. It was an 
impressive unexpected result to find that the ratio of 
viscosity to entropy density in the Af — A SYM plasma 
at strong coupling turned out to be an almost universal 
prediction, whose value comes closest to the results that 
fit the data obtained from RHIC and LHC. 

The aforementioned breakthrough boosted an enor- 
mous activity, and opened up a very interesting field 
of research that attracted both physicists from string 
theory as well as from heavy ion and QCD communities. 
The core of this program is to obtain the physical 
transport properties of a relativistic quark gluon plasma 
(much of the research in this field is discussed at length 
in the recent review in Ref. 3] which also contains an ex- 
tensive reference list). In its more crude form, the initial 
attempts were really only dealing with adjoint degrees 
of freedom. To say that A'c = oo is a sensible approxi- 
mation to Nc = 3 is probably not as severe a statement 
as to say that in real QCD the number of flavors is zero, 
or even negligible as compared with the number of colors. 

From the holographic perspective, adding flavor 
degrees of freedom that transform in the fundamen- 
tal representation of the gauge group is fairly well 
understood in the so called quenched approximation. 
A quenched quark, a term borrowed from the lattice 
literature, means a non-dynamical quark. In pertur- 
bation theory diagrams with internal quark loops in 
it are neglected. The easiest realization of this is the 
case of very massive quarks. For light quarks there 
is another quenching mechanism in the so called 't 
Hooft limit,E' Nc ^ oo and A = gyM-^c flnite. In 
such a limit, diagrams with i nte rnal quark loops are 
suppressed by factors of iVj/iVc.'^ Therefore keeping Nf 
fixed induces a quenching effect even for massless quarks. 

From the gravity dual perspective, this approach 
amounts to having a finite number of probe flavor branes 
which do not backreact onto the geometry that the color 
branes create. There have been several attempts to 
go beyond the quenched approximation by considering 
metrics sourced by an infinite number of Nc color and 
Nf flavor branes. Apart from the works explicitly 
discussed in the rest of this paper, we refer the reader 
to Refs. [7] and [HI and to the excellent review in Ref. 
[HI for an extensive list of references on the topic. The 



present review will be concerned exclusively with one 
such strategies, which we will term as "smeared flavors^'' . 
Other proposals where the backreaction of a localized 
stack of D7-branes is taken into account are included in 
Ref. \M 

Outline: Our aim in this review is to provide a very 
general account of the so called D3-D7 holographic quark 
gluon plasmas, includ ing fin ite temperature and chemical 
potential dependence. f^^^^l^ In section II we review the es- 
sential ingredients of the "smeared flavor" construction, 
in the particular setup of D7-branes in the background of 
D3-branes. The geometry sourced by the D3-branes will 
be of the general type having a Sasaki-Einstein internal 
space. In section III we introduce the ansatz that incor- 
porates the desired physics, and the equations of motion 
are derived. Section IV is devoted to solutions to these 
equations. An analytic solution is possible as a pertur- 
bative expansion in small values of e ~ Nf/Nc and the 
baryon density Uf,. In essence we extend and fill in some 
gaps of the work presented in Ref. [14] and provide the 
explicit solution up to second order in e and first nontriv- 
ial order in rif,. In section V we shall show how to obtain 
consistent thermodynamics from this solution and dis- 
cuss the susceptibilities. Finally, in section VI we explore 
some other physical predictions of the mentioned model. 
In particular, the hydrodynamic and optical properties 
of the quark gluon plasmas are examined. We conclude 
with some appendices containing complementary mate- 
rial. 



II. D3-D7-BRANE SYSTEMS 

Introducing a stack of Nc D3-branes at the tip of a 6d 
Calabi-Yau cone over a 5d Sasaki-Einstein manifold 
provides a gravity theory dual to a supersymmetric field 
theory in 4d which preserves (at least) Af = 1 supersym- 
metry (see Refs. [151 and [TBI for a review see Ref. [T7|) . 
By far the most studied example of this setting is when 
the Sasaki- Einstein manifold is S^, so that the gravity 
theory is that of AdSz x S^, and the dual field theory is 
the conformal field theory of A/" = 4 SYM with SU{Nc) 
gauge symmetry. However, these theories do not yet in- 
clude matter in the fundamental representation of the 
gauge group; to do this, we must deform the theory by 
adding the fiavor degrees of freedom. This section re- 
views briefiy the probe approximation and the smearing 
approximation, two ways in which fiavor can be added to 
the field theory. 



A. Introducing Flavors 

We can add matter in the fundamental representation 
of the SU{Nc) gauge group by introducing D7-branes 
into the system; a stack of Nf coincident D7-branes 
introduces a U{Nf) global fiavor symmetry into the 
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theory. Then, we have three types of open string 
excitations: open strings ending on two D3-branes which 
describe gluons in the adjoint of SU{Nc), open strings 
ending on both a D3 and a D7-brane which describe 
quarks in the fundamental of SU (Nc) and open strings 
ending on two D7-branes. The latter are in the adjoint 
representation of U{Nf) and thus naturally represent 
mesons (for a review, see Ref. ITS)) . 

Adding D7-branes in this fashion was first done for 
TV = 4 SU{Nc) SYM in Refs. [El EQl EH The introduc- 
tion of flavors in this way breaks the Af — 4: supersym- 
metry to TV = 2 and also conformality in the field theory. 
The following intersection of the D3 and D7-branes is 
used 
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Here, a^~^ are the coordinates on the internal compact 
manifold . The D3-branes are localized along the AdS 
radial variable r at r = 0. If we want massless quarks, 
the D7-branes should extend to the origin and intersect 
the D3-branes; for massive quarks, the D7-branes should 
have a radial profile and only extend down to a certain 
radial position - this position is then related to the 
mass of the quarks. We will be mainly concerned with 
massless quarks in the following. 

The contribution to the full action of the flavor D7- 
branes compared to the contribution of the D3-branes 
can be calculated a^ 

^D7 ^ Cd3 , (1) 

where Cd7,D3 are the respective Lagrangians and we have 
used the 't Hooft coupling A = gyM^c- Essentially, this 
implies that the parameter Nf/Nf, (or, to be more pre- 
cise, XNf/Nc) is the measure of the backreaction of the 
D7-branes on the AdS background geometry sourced by 
the D3-branes. In fact, the backreaction of the Nf D7- 
brane s is ze ro in the 't Hooft scaling limit (keeping Nf 
fixed) ™l 



i.e. when we take the number of flavors to be very small 
compared to the number of colors. This is called the 
probe approximation for the D7-branes, as they probe 
the geometry sourced by the D3-branes without deform- 
ing it. From the field theory perspective, we are ignoring 
the quantum effects caused by the quarks - this can be 
seen by considering different diagrams contributing to 
a given physical process: the contribution of a diagram 
with k internal quark loops compared to a similar 
diagram with none scales as [Nf/Nc]^. Thus, using the 
probe approximation for the D7-branes gives us quarks 



in the dual field theory which are essentially external 
non-dynamical objects which do not run in loops. This 
is the quenched approximation for the fundamental 
matter in the field theory, a term borrowed from lattice 
literature; the approximation becomes exact in the limit 
§. 

The dynamics of the D7-brane in the probe approxi- 
mation is governed by its DBI action, 

Sdt^-Tj J d«Ce-*^det(G,b + B^b + ^^a'P^^) , (3) 

where Gab is just the puUback of the AdS^ x metric 
sourced by the D3-branes. In the probe approximation, 
one only needs to consider this action to determine the 
full dynamics of the D7-branes. 

B. Finite Baryons Density 

QCD is expected to show a rich phase diagram at finite 
baryon density, including a color-flavor locking phase 
at high baryon density and relatively low temperature. 
In addition, neutron stars are astronomical objects ex- 
pected to be at virtually zero temperature but non-zero 
baryon density. Thus, it would be interesting to use 
holography to probe properties of strongly interacting 
field theories like QCD at finite baryon density. 

In QCD, the subgroup U{1)b of the U{Nf = 3) global 
flavor symmetry is identified with the baryon number. 
Thus, we wish to introduce a finite U{1)b density (J*), 
where is the U{\)b current. This baryonic current is 
dual to a vector field A^ living on the D7-branes.^^^^ 
The precise correspondence is, as calculated in the probe 
analysis!^ 

At{r) ^Oq = + +1 {q^Vtq - {Vtq)U) 

+l{q{Vtq)^ -Vtqq^), (4) 

where iPt'4' are the fermions in the (anti)fundamental 
representation of SU{Nc), and q,q are scalars in the 
(anti)fundamental representation; V is the covariant 
derivative in the correct representation; Og is the quark 
number operator. Note that we only need the time 
component of the gauge field At , and that this field only 
depends on the AdS radius r. A complete analysis of the 
stability of the system in the probe approximation can 
be found in Ref. |25| while recent studies for imaginary 
chemical potential are reported in Ref. 

In the probe analysis, the asymptotic behavior of At{r) 
for large r iJ^ 

Mr)-^^-^, (5) 

where is the quark chemical potential and a is propor- 
tional to the quark density Uq. The quark density Uq 
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is related to the baryon density Ub as ub = riq/Nc as 
there are Nc quarks in a baryon; analogously /is = fJ-Nc- 

In Ref. 27, the so-called linearized backreaction on the 
supergravity geometry of including a non-zero At{r) is 
calculated. The result is that a component of the RR- 
field i^(3) is turned on (see also Ref. For the zero 

temperature case, it is found that F(3) = Fi23dx^ A dx^ A 
dx^, where 



23 



Nr 



(6) 



Here c? is a constant related to the quark density as 
d = nq/Vo\{S'^) and is obtained as a constant of motion 
from SCm/SFrt, with Frt = 5^ A* .1221221 This constant 
f 123 component of i^(3) will be an important motivation 
for the choice of ansatz for F(-3') below. 



C. Beyond Quenched Flavors 

The discussion in this section is mainly inspired by 
that in Refs. [5] and [2^ We shall keep the discussion 
fairly minimal here, we refer the reader to Ref. [9] 
for an in-depth review of the smearing procedure for 
backreacting D7-branes. 

As discussed above, in the 't Hooft scaling limit ([2|, the 
backreaction of the D7-branes can be neglected. How- 
ever, if we want to introduce dynamical flavors that run in 
loops in our field theory, this scaling limit is insufficient; 
we need to take the backreaction of the flavor branes into 
account. We can do this by using the Veneziano limit,-^ 



Nf 



Nc 



const. 



(7) 



Thus, now we keep the ratio Nf/Nc fixed while taking 
both Nf and N^ to infinity. The ratio Nf/Nc essentially 
determines the magnitude of the backreaction of the 
flavor branes, as we will see. 

Clearly, we can not neglect the backreaction of the 
Nf branes on the AdS geometry in this scaling limit. 
Now, instead of considering the AdS background as a 
solution to the HB supergravity field equations and then 
minimizing the D7-brane DBI action in this geometry, 
we will minimize the full supergravity action Ssugra 
together with the D7-brane DBI and WZ action Sdj 



S — Ss. 



UGRA 
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D7 ■ 



(8) 



The D7-brane action Sd7 is an integral over the 8d world- 
volume of the D7-branes. Such integrals are typically 
very difficult to deal with; this is why the powerful tech- 
nique of smearing was introduced, which we will now 
discuss. 



D. Smearing 

Let us consider the Abelian DBI action (without form 
fields) , 



(9) 



Sdbi= / rf«eV-G, 

J Ms. 



which is an integral over the 8d world volume of the D7- 
brane. To write this as an integral over the entire lOd 
space-time, we can use a localized two-form VL2 such that 



-a A 13, 



(10) 



where a, (3 are one-forms which are orthogonal to A^s f^nd 
g refers to the ten-dimensional metric. Then we have 
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Note that |f^2| = EJ^^'-*! (where f^^'^ are the decom- 
posable parts of f22) is the coordinate-invariant modulus 
of 0.2 , which in the supersymmetric or trivial embedding 
case reads 



0| 



1 



21 "a/a^^p^ 
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(12) 



Since the (5-functions in ( 10 1 are virtually impossible to 



work with in the general case, it would be advantageous 
to get rid of them; this is precisely the purpose of the 
smearing procedure. One can effectively view this as 
taking a large amount of D7-branes and smearing them 
over their transverse space.f^Sl 

An apt analogy for this smearing procedure using 
electrically charged lines is given in Ref. |9l Let us here 
give our own short analogy. Consider a disc containing 
a number of points. If the number of points in the disc 
is small (see fig. [T]), then we safely approximate the disc 
as having nothing in it - this is comparable to the probe 
approximation, where the number of D7-branes (the 
points) is small compared to the number of D3-branes 
(the disc's area). However, when the number of points is 
large (see fig. |2| , then we can no longer purport the disc 
to be "almost" empty - we have to take the points into 
account in the description of the disc. Taking each indi- 
vidual point into account would be quite a chore, indeed. 
However, if the points are relatively homogeneously 
spread throughout the disc, then we can approximate 
the distribution of points using a continuous distribution 
function - this turns our complicated disc in fig. [2] into 
the smooth, more manageable, disc in fig. [3] We have 
now effectively "smeared" the points in the disc. 

It is in any case obvious that such smearing procedures 
are not by any means new in physics. For example, a 
liquid is really a collection of small particles, but from 



5 





FIG. 1. A disc FIG. 2. A disc FIG. 3. The large 



with a small num- 
ber of points in it. 



with a large num- 
ber of points in it. 



number of points 
in the disc have 
been 'smeared' 
out into a uniform 
distribution. 



a macroscopic viewpoint we can 'smear' the particles 
and just consider the fluid as a continuum. Note that 
for the fluid, it is important that there are a large 
number of particles; the situation is the same when we 
consider D7-brane smearing - we need a large number 
of them before the smearing procedure gives us a good 
approximation. 

The two-form as deflned above essentially tells us 
where the D7-branes are located. Without smearing, it 
consists of delta functions and is not easy to work with; 
the smearing procedure tells us how to 'smear' this two- 
form. The crucial observation in Refs. [311 IMl is now 
that a smeared ^2 is not arbitrary if we want to preserve 
supersymmetry. It turns out we require, for massless 
flavors, 



Q2 — —'^QfJRE : 



(13) 



where Jk e is the Kahler two- form of the Kahler- Einstein 
base of the compact Sasaki-Einstein manifold X5, and Qf 
is a constant related to Nf as 



Qf 



Vol(X3)5,7V_ 



/ 



4Vol(X5) 



(14) 



where Vol(X5) is the volume of the internal manifold, 
and ¥01(^3) is the volume of the submanifold X3 of ATs 
that the D7-branes wrap. Note that every D7-brane in 
principle has a different manifold which it wraps, but 
all of these Xa's are related to each other by symmetry 
transformations, so Vol(Ar3) and thus Q/ is well-defined. 

Practically, we use the smeared Q2 by making the fol- 
lowing substitution for any eight-form X(^gj which we in- 
tegrate over the D7-brane worldvolume 



X 



Ms 



(8) 



Mio 



X(g) A ^2 



(15) 



In particular, for pull-back quantities such as in the Wess- 
Zumino part of the D7-brane action, we have 



C, 



(8) 



C(8) A Q2 



(16) 



We note that can be seen as a magnetic source for 



the form-field Fi 



(1) 



dF^i) = -n2 



(17) 



since the D7-branes source an eight-form potential. This 
means that we can take 



= Qf{dT + Ake) , 



(18) 



where r is the U{1) fibration coordinate of the Sasaki- 
Einstein manifold X5, and Ake is the Kahler potential 
{(IAke = ^Jke)-^^ We also note that we can give the 
quarks mass by giving the D7-branes a non-trivial profile 
in the radial direction; we can do this by introducing a 
function p{r): 



F^,)=Qfp{r){dT + AKE) 



(19) 



Of course, p{r) — > 1 for r 00, as the quarks are effec- 
tively massless in the UV at energies much higher than 
the quark masses; and p{r) — for r < Vq for a certain 
radius Vq as the quarks can be effectively integrated out 
at energies much lower than the quark masses. The func- 
tion p{r) can be calculated explicitly and can be shown 
to be a smoothed-out theta-function.l^ 



E. Validity of the smearing 

The smearing technique described above makes use 
of the Abelian DBI action in order to take into account 
fiavor effects. This approach has some limitations; we 
discuss in this section three topics which sharpen its 
regime of validity. 



As a first point ,1^21231 the DBI action has been shown 
in Ref. f^Slto take into account the leading contribution 
of the open string sector, i.e. the 0{gsNf) terms in the 
topological expansion of the string amplitudes. This al- 
lows to resum, in the dual field theory, the leading effects 
in Nf/Nc, the "one- window" graphs in the language of 
Ref. 5, i.e. the graphs containing one loop of matter in 
the fundamental representation ("quark loops"). The 
result is surely an improvement with respect to the probe 
limit, which does not include any loop of fundamental 
matter. It does not provide the full backreaction of 
fiavors, though. The latter would require to resum 
the contributions from the graphs with any number of 
quark loops. Going beyond the "one-window" graph 
contributions is surely an interesting issue; a proposal 
can be found for example in Ref. fMl Nevertheless, 
if gsNf is small, the smearing technique provides the 
leading fiavor contributions, since the graphs with k 
quark loops are weighted by powers of (gsNj)''. As we 
will see in the following, we will indeed require gsNf to 
be small. Note that in any case the effects taken into 
account by the DBI are quite non-trivial, non-linear 
in gsNf because the "one-window" graphs are resummed. 
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On top of this, the use of the Abelian DBI is fully 
justified only if the minimal distance between two 
branes is much larger than the string length, i.e. the 
strings connecting different branes are massive and 
can be ignored among the light modes. If R is the 
typical size of the n internal directions along which the 
smearing is performed, the minimal distance between 
two branes will be of order R^/Nf. This poses the 
restriction Nj ^ _R". In the example of this paper, 
^ \,n — 2, and considering also the requirement 
N f ^ 1, it turns out that the smearing approximation is 
justified in the parameter window 1 ^ Nf <C Since 
in this window the non-diagonal modes among the D- 
brane spectrum are essentially decoupled, the smearing 
technique describes the flavor-singlet sector of the dual 
field theory. An improvement in this direction could 
only be captured by a consistent non- Abelian DBI action. 



To summarize, smearing the D7-brane action is a pow- 
erful technique which simplifies the D7-brane action to 
only depend on the AdS radial coordinate; our equa- 
tions of motion will (usually) thus be ODEs instead of 
PDEs, and in any case result in equations without 5- 
functions. However, in addition to the validity limita- 
tions as described above, we must also beware that the 
smearing is only a good approximation if the number of 
flavor branes is large and moreover they are uniformly 
distributed throughout the directions transverse to the 
branes. So, the global U{Nf) flavor symmetry sourced 
by Nf coincident D7-branes is broken to U{l)^f when 
we smear them, as smearing requires the D7-branes to 
sit at different points in the transverse space. 



III. SET UP & ANSATZ 



Finally, let us point out a limitation concerning 
the use of the DBI in the flnite temperature case. It 
has been known for a while (for recent work in this 
direction see for example Ref. that the DBI does not 
capture the thermalization effects of the theory on the 
world-volume of the D-branes: its energy momentum 
tensor does not coincide with the black brane one. As 
such, it constitutes only an approximation to the dual 
description of a fully thermalized gauge theory with fun- 
damentals. The effects that are not taken into account 
scale as Nf/Nc times some power of the temperature (in 
the zero temperature limit the discrepancy is of course 
absent). So, the DBI is a good approximation in the 
probe approximation or at small temperatures. Thus, 
we will be describing the small temperature regime of 
the dual field theory; this limitation is anyway imposed 
upon us also by another reason: the presence of a Lan- 
dau pole in the UV in the theories we will be considering. 



We wish to study a field theory which is dual to the 
Type IIB gravity theory on AdS^ x X^. This field theory 
must be deformed by Nf dynamical fiavors, so we will 
add a large number of D7-branes to the gravity theory, 
smeared in the way discussed above. Finally, we turn 
on a finite temperature and finite baryon density in this 
field theory; in the gravity theory this translates as the 
introduction of a black hole in the geometry and a gauge 
field At on the D7-branes, respectively. 

The method we must use to study this gravity theory 
is as in Refs. [TT] and [TH i.e. the second order equations 
of motion must be solved for. We can not use supersym- 
metry considerations as we are interested in the system 
at finite temperature so that supersymmetry is broken. 

The starting point is thus the full supergravity action 
together with the D7-brane action. This action is 



S 



D7 



(3) 



^2* p2 



C(4) A i/(3) A i^(3) -I- Sfi , 

d^Xe^\/- dct(G + e-*/2j-) + f C^e'^ . 



$ 1^3 



1 



(3) - 4^(5) 



(20) 
(21) 



Here we have defined J- :— 27ra'F + B. Note that 
K) = (also for i?(3)). The 

form-fields F(^p^ are the gauge invariant field strengths of 
the potentials C(p_i) (eventually shifted by -ff(3)). Note 
that the potentials C(p) are not well-defined in the pres- 
ence of electric and magnetic sources on the D-branes. 
The self-duality of -^(5), i.e. i^(5) = *F{5), is not deriv- 
able from this action but needs to be imposed by hand. 
The equations of motion that extremize this action are 



Einstein's equations for the metric, the equation of mo- 
tion for the dilaton, and the equations of motion for the 
form-fields (which are derived and stated in full in Ref. 



A. Ansatze 

For the metric, we concentrate on the convenient 
ansatz 



ds^ = h{a)-^/^[-bia)df + dxj] + h{ay/^ [b{a)Sia)^F{ayda^ + Siafds^^ + F{af{dT + Ake? 



(22) 



Our radial variable is a (with dimensions of length^). 
The ansatz above is suitable for a charged black hole in 
a deformed AdS space times a squashed Sasaki-Einstein 
manifold. The squashing is produced by the backreac- 
tion of the D7-branes trivially embedded in the geometry 
(i.e. corresponding to massless flavors). 

Having the metric functions b, h, S, F depend only 
on a radial variable a is possible due to the smearing 
procedure outlined above. In addition, we also have 
the dilaton $(cr) which we also take to depend only on a. 

As we saw above, to introduce a finite baryon density 
in the dual field theory, we use the following ansatz for 
the gauge field on the D7-branes: 



A = At{a)dt. 



In addition, we set B 



(3) 



(23) 

throughout. This 



also implies that T = 2-Ka'A[da A dt and T A J- = 0. 
The equations of motion and Bianchi identities which 
determine the Ramond-Ramond form-fields then become 



d(e2**F(i)) 
d(e**F(3)) 

die-" * i/(3)) 



dFt^s) = 



0, 
0, 

= e*F(i) A *F(3) 



Fi5) ^ F(^3) 



(24) 
(25) 
(26) 
(27) 



det(G + e-*/2j-),5(2)(i?7), 

(28) 
(29) 



Note that in (27 1, the final term is not explicitly written 
as a form, following Ref. [351 Finally, the equation of 
motion for the worldvolume gauge field At coincides 
precisely with the derivative with respect to a of the 
right hand side of (27 1. Hence with our ansatz having 
B — the equation of motion for At is implied by the 
equation of motion for -ff (a) • 



To solve these equations of motion and Bianchi iden- 
tities, it suffices for F^i) and F(5) to assume the ansatze: 



■^(1) 



F 



(5) 



QfidT 
Qc{l + 



f Ake) , 
*)e(X5) ., 



(30) 
(31) 



where £(^5) is the volume form of X5 and we have in- 
troduced the constant Qc through 



Vol(X5 

(2^)4g,c 



(32) 



Note that at finite temperature, we no longer have 
supersymmetry to motivate the D7-brane embedding as 
in (11); however, it is checked in Refs. [TTJ[T3]that the 



trivial embedding is still consistent for the D7-branes for 
massless flavors. We will only consider massless flavors 
in what follows, see Ref. [TT|for an example of equations 
concerning massive flavors at finite temperature. 

Finally, we are left with the form- field F(3) . Motivated 
by the discussion above, we expect that introducing a 
non-zero field At on the brane will imply that F(3) in- 
cludes a non-zero component ^123^2;^ A dx^ A dx^, where 
we postulate F123 to be a constant. However, this is 
not the end of the story; an analysis of the equations of 
motion for the form-fields shows us that consistency re- 
quires the introduction of additional components of i^(3) . 
We can motivate our choice by noticing that there is a 



WZ term Cjg) A in the action (|21|); as soon as Cfg) and 



^(6) 



are non-zero, this term must be taken into account. 
Let us take the following ansatz for the C(6) potential: 



C(6) = J(cr)dx'"^ A (dr + Ake) A . 



(33) 



This term gives us two extra terms in our ansatz for 
F(3) (through F(7) = dC(6) and F, 



^(3) 



* F, 



{7)) 



Our full ansatz for _F'(3) then contains three terms in total 



i^(3) = Fi23dx^ A dx'^ A dx^ - dtAVL2+ SQfJe^'^bF'^ dt A da A [dr + Ake) 



(34) 
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B. Equations of Motion 



Once we have specified our ansatz, we are left witli 
the metric function b{a), h{a), F{(t), S{a) as well as 
the form- field function J {a), the gauge field At{a) 
and the dilaton $((t); these will be uniquely deter- 
mined by the equations of motion and the constant -^123, 
which is proportional to the baryon density as we will see. 



First of all, At{a) and J(cr) are determined in terms of 
the metric functions by the two equations (27 1 and (29 1, 



which give us 
2Tr A[{a) = 



{QcFi23 + 8Q}J)bS''F 



WQjS^F'^ 



-8 J e^'^bF' 



e-'^{QcFi23 + SQ}jy 
-<s>jn 



d 
da 



(35) 



Then, we need the equations of motion for the met- 
ric functions and the dilaton. We can obtain them by 
calculating the left- and right-hand sides of the Einstein 
equations and equating them, or by directly inserting our 
ansatze into the action (20), (21). In the latter case, we 
are left with an effective one-dimensional actionl^l 



yol{X5)V, 



1,3 



2k^q 



Lidda , 



(36) 



where Vi^3 denotes the (infinite) volume integral over 
the Minkowski coordinates, and 



nd = -l (log' h)^ + 12(log' Sf + 8 log' F log' 5 - ^$'2 + !^ (log' h + 8 log' S + 2 log' F) - 4Q} ^ 

~ ^^^^^^ + 2AbF^S'^ ~ ^Fi^^e'^bh^F^S^ ~ ^Qje'^^bS'' - Ae^'^FQfS''^ -{2T:a'A'^Y + eH^F'^S^ 
-i2Q)e-'^hF^J^ ~8Q){2na'A't)J . (37) 



Since At only enters the Lagrangian with its first deriva- 
tive, it must be connected to a constant of motion. In 
fact, this constant must be set as follows: 



5A't 



2na'QcFi 



23 : 



(38) 



relevant equations of motion and obviously coincide with 
those obtained directly from the Einstein equations.^'* Fi- 
nally, the equations of motion as derived from ( 39 ) must 



be supplemented by the so-called zero-energy condition 
iJ = 0, with 



as this gives us exactly equation (35). In fact, we can 
use (35) to eliminate A'^. in favor of F123, obtaining the 



Lagrangian 



Lid — L 



Id 



6Lid 



(39) 



The Euler-Lagrange equations from ( 39 ) then give us the 



H = -Lid + J2^^%^^ ^^ = {b,h,F,S,'^,At,J}, 

i * 

(40) 

which is a first-order differential equation in the un- 
known functions. 



In any case, the end result is the following set of equations of motion (also restating those determining At and J 
from above) 
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{log by = be^F^F^^^h^S'' 
(log/i)": 



4(27raM;)2e*/2FQ/S'2 



^254 



2 y3(wl02TWF258 



4e-*g2j'2 
^25"* 



326e-*i^2Q2j2 



(log 5)" 
(logi^)" 



2^ a; (a) 



52g3*/2^3Q^^l 



^-(27ra'A',)2 + 625*^25-8 



166e-*F2g2j2 



(41) 
(42) 
(43) 



4c4 i,„'S> ip2 ip2 u2as ^be'^'^Q'jS^ (27ra'Ai)2e*/2^Q^5'2 



1, 



ibF^S^ - -be^F^Ft^^h^S- 



16&e-*F2Q2j2 



v/-(27ra'yl02 + 626*^^25-8 



2e-*Q2j/2 
F254 



(44) 



;6e*F2^^2^3/i2 5« + 6e2*g2.58 + 
-326e-^2Q2^2_4e^:/- 

(Qci^l23 + 8QU)bS^F 



IGQjS^F^ + e-^{QcFi23 + 8Q}jy 
{QcFi23+8Q)J)bS^F 



'{QcFl23+8Q}j)^ 



2fc2g3$/2j.3g^5l0 

v/-(27ra'A^)2 + 62ei'i^2 5.8 



+ 8J(e-*6iJ^2 



16Q254^2 



together with the first-order constraint 



2e*/2FQ/5Y-(2^a'^i)2 + 626*^^2^8 

(45) 
(46) 

(47) 

1 



= log' /i log' b + ^(log' hf - 12 (log' Sf - 4 log' 6 log' 5* - log' 6 log' F - Slog' log' S + ^^'^ 
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- ibF^S'' + 2AbF^S^ - 36*6/^2^^258 _ l;,e2*Q2 58 _ ^^'^'^^'^'^^'^ 

2 123 2 ^-^ ^-(27ra'A02 + 62e*F258 



4e-*g2j'2 
^254 



-326e-*i^2Q2j2 



(48) 



Now that we have obtained the equations of motion, we 
can finally proceed in order to find their solutions. 



IV. THE PERTURBATIVE SOLUTION 

In order to find a complete solution of the equations 
of motion, we should in principle proceed by numerical 
integration. Fixing the asymptotic behaviors one should 
glue them numerically so as to determine the integration 
constants. This kind of analysis is especially difficult in 
the case at hand due to the presence of a Landau pole, 
and stability of the numerical integration becomes a 
very relevant issue. 

The appearance of a Landau pole is expected on 
physical grounds as we are perturbing conformal field 
theories, as e.g. N — A supersymmetric Yang-Mills, with 
the addition of fundamental matter. The one-loop beta 
function will be positive and proportional to Nf and 



the running coupling constant will have a divergence at 
some UV scale Kuv In the holographic dual it will be 
seen as a divergence of the dilaton field. As is the case 
for QED, the meaning of this is that we are dealing with 
an IR effective theory, and a suitable completion should 
replace it in the UV before the Landau pole is reached. 



To gain information about the full solution, a clever 
approach is to follow a procedure analogous to the one 
considered in Ref. |37]to construct an approximate black 
hole solution of the Klebanov-Tseytlin model. Given 
the fact that, for M the number of fractional branes 
and N that of regular D3-branes, the exact black hole is 
known when M = 0, Ref. ISTl considered a perturbative 
expansion with parameter IvP/N ^ 1. Already at 
first order, such perturbative analysis gave notable re- 
sults. The same strategy, supplemented by a numerical 
analysis, has been applied in Ref. 1381 where a solution 
at finite charge density for di-baryons has been discussed. 
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In the flavored setup, the natural parameter for a 
perturbative expansion is given by the effective couphng 
gprpNf which weighs the vacuum polarization effects due 
to the dynamical flavors. In the unquenched smeared 
D3-D7 model, the supcrgravity plus DBI solution is reli- 
able in the Veneziano limit Nf ,Nc^oo with XNf /Nc = 
constant, provided the 't Hooft coupling A = g^rpN^ is 
large. This in turn implies that Nf/Nc <^ 1 and the 
effective coupling gpT^f ~ XNf/Nc can be of order one. 
Inspired by Ref. [37] we will now consider the perturba- 
tive regime in the Veneziano parameter e ~ XNf /Nc 1. 

In this section we will derive an analytical perturbative 
solution at finite temperature which has a nontrivial 
backreaction of the flavor degrees of freedom, as well 
as a nonzero baryon chemical potential. The price of 
having an analytic solution is that we must keep the 
deformation small, so that we can perform a perturbative 
expansion around the original unflavored theory. 

In a first approach we expect roughly three indepen- 
dent parameters, which we will name rh,e and S. The 
first one is related to the temperature, r/j ^ T. For 
e = S = it will parametrize the usual AdS^ black hole 
solution with T = r^/ (jry/Xa'). The second one, £ ^ Qj, 
is a measure of the backreaction of the flavor degrees 
of freedom and, in fact, turns out to be related to the 
flavor-loop counting parameter. Finally, related to the 
baryon density we will have another parameter d ^ A'^ . 



Concerning e, it turns out that the relevant flavor 



counting parameter is scale dependent. The best way 
to think of it is to notice that, in the equations of motion 
^7] ), at zero baryon density (nf, = J = = 
= 0), Qf and the dilaton $ always appear in the 



(41) 



^123 



combination Q/e*. Hence, consider a natural split of the 
dilaton fleld at an arbitrary "anchoring" scale A* ^ cr*, 
$((7) — +<P{(^) where 0(ct*) = 0. Then at that scale 
is deflned by 



Qr 



Vo\{X3 



167rVol(X5 



-A. 



Nj_ 

Nr. 



(49) 



where A* — ingge " Nc and we have made use of (14 1 



Hence weighs the backreaction of the D7-branes and, 
in fact, can be read as a flavor- loop counting parameter 
in the dual field theory. Keeping small is tantamount 
to imposing a large separation of the scale and the 
dangerous Landau pole aLP^ 

In order to derive an analytical perturbative solution 
we will take both e* and d to be much smaller than one. 
All the functions in the ansatz will be expanded in power 
series of e™(5". To get started let us introduce the follow- 
ing leading scaling behaviors 



23 



Qf 



(50) 



The reason behind this choice goes as follows: inserting 
these expressions into (41 )-(47l and recalling the dilaton 
splitting as $(a) = ^(cr), with 0(cr*) — 0, one 

readily arrives at the following system of equations 



(log b)" = iej^ y+6i el5^ e"^ hF^P + 8 el5^ e 



(logft)" 
(log 5)" 



5-4^2 



X 
Y 



2 ri6 



J' 



f2 



F^S^ 



32 eiS'' e-f bF'j'^ + 4 elS^ e^"^ 



J'2 



^2 5*4 



b'^F^S 

Y 
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16e^52e-0 5^2j2_^2j2 



(logF)" = ^bF^S^ - \ele^t' bS^ -e,5^^ + lQ el5^ e"^ bF^P - 2 el5^ e"^ 



(51) 
(52) 
(53) 
(54) 



2e*e"*/2FS'2y - 326^52^-0^^2 j2 _ Ae^S"^ €'">' + e*<5^ ' ^^^^ 



(H-8e*J)6FS'4 



-Hf'^j . 



(56) 



The constraint equation ( 48 1 reads 



= - ^ log' h log' b + i (log' hf - 12(log' Sf - 4 log' b log' S - log' b log' F - 8 log' F log' S +]^(t)''^ 



-^-^^AbF^S^ + 2AbF^S'' 



4,30/2^2 j.3gl0 1 

- Y *2 



eiS' -32be-fF^J^ 



(57) 



4e-^ J'2 1 ^ 
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In these equations X, and Z are shorthand for 



X = 



(l + 8£,J)^e^/"5^f3gio 
16F254 + (52e-'^(l + 8e*J)2 



(58) 



(52(l + 8e,J)262f2g8 



Y = X b'^e't'F^S^ - 



r(59) 



(60) 



however, in order to stress that our parameter is going 
to be perturbatively small. 

The strategy now is to set up a perturbative ansatz in 
both e* and i5, by means of which, all functions admit an 
expansion like 



i,j>0 



(63) 



The system ( 51 1-( 56 ) allows for a systematic expansion 
of all the functions in powers series of e* and S^. This 



is essentially the main effect of the scaling relations ( 50 1 



Once all the functions have been solved for, the world- 
volume gauge field can be obtained from the following 
relation 



{l + 8eJ)bFS* 



16F2S'4 + 52e-0(i + 8e*J)2 



(61) 



which is already first order in S. From this expression it 
is fairly obvious that J(cr) is a higher order effect in 
and hence arises as backreaction of the baryon density 
onto the flavor branes. 

Actually, it is convenient to introduce a dimensionless 
parameter S defined as 



4,3 



(62) 



where the factor 4 is introduced in order to make it 
precisely 6 = d oi Rcf. 23, We keep the Greek symbol. 



The reader may have noticed that we have changed 
radial coordinate cr — r. This is done after integrating 
the equations of motion in a and the new radial coor- 
dinate "gauge" is chosen to make the dx^ component 
of the metric equal to r^/R^, just the same as in the 
unflavored solution. This allows for an easy comparison 
of the two situations. 



Equations (51)-(56) are solved for order by order in 



£*, 6. At each order, two conditions must be imposed in 
order to fix the two integration constants that appear 
in each function. One of them will be regularity at 
the horizon r = r^- The other one will be matching 
the supersymmetric solution found in Ref. [29] at some 
UV energy scale r = r^. In summary we have three 
energy scales: the horizon radius r^, the dilaton an- 
choring point r^, and the UV cutoff scale rg. and 
Ts are arbitrary, and in Ref. |TT]they were fixed to — r^. 

A solution containing the three parameters is un- 
wieldy, and we sketch here just an example up to ©(e*) 



0(r) 



1 + 6 log — 



2r4 



.4^ 



2rt - 4 



(64) 



which naturally satisfies 4>{r^,) = 0. In Ref. [TT]a solution 
to order was presented in equation (2.22) with 
— Vs- Namely, the scale where the 't Hooft coupling is 
defined and the one where the thermal solution matches 
the supersymmetric one coincide. The solution is well 
defined in the range < r < r*. If one is interested 
in IR quantities, like thermodynamic properties, or 
transport coefficients, it is more sensible to work with 
a solution where all the running parameters are defined 



at the relevant scale, e.g. with eh = Q/e*'- The 
solution with = and Ts — > oo was used in section 
3 of Ref. [n as well as in eq. (2.7) of Ref. [TH In Ref. 
[13] we presented an extended version of the solution 
which contained corrections up to order S^. In that 
paper, rather than the metric coefficients, we gave the 
expression for the functions in the ansatz, F, S, etc. 
We will provide here a final form for the metric, dila- 
ton and At potential. Namely, writing the line element as 



ds^ = Gttdt^ + Gxxdx\ 



G^rdv 



GrridT + A^ 



(65) 
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we find 

Gttir) 

Gxx{r) 
Grr{r) 



\{2r^ - rl) + rV4(log 2 - 1) + r2(r^ - 2r^) log (l + ^) 
2r4R2 



B? 1 

^4 



-34r*r^ + 5r^ + rV^(29 - 171og2) + 17r2(2r* - r^) log (l + ^) 

24r4R2 



+ 



(66) 
(67) 



■4 



-5^ 
+..., 



8.2 J. ^ 7.4^6 _ .10 j_ .2.8 (^^g 4 _ 3) _^ ^2^5^.8 _ 9^.4^4 ^ ^8) (l + ^j' 



4+*^ ( 4r4r2(r4-r4) 



11 -24 log 
96 



-42r^r^ + 21r*5r^ + 49r^r^ + 5rj^° + (-33 + 14 log 2) + Ir^iQr^ - Qr^r^ + r^) log (l 



24r4(r4 - r^) 



1 ,5-241og- , I , 



(68) 



2.4 



-tIt'G{t) + ( 2rV^ - 3r^ + 3(-2r^ + r^) log ( 1 + ^ 



G..(r)=i?2<;i-e,^+4' 



3 + 8 log 



96 



,52 ('^ + 34-, ^ 
20 40 



(69) 



2.4 



-2ryG{T) + (22r^ - 3rV^ - 5r^) + 3r2(-2r^ + r^) log 1 



(70) 



(/>(r) = e?, log — + £/, 



^ (l + 31og— ) log( — 
6 V "Th \rh 



' ?r4I-2,-^-i¥+OW-291„g2+('lI + l»r! 



120 



log fl 



1 + ^1+ J^Li2 ( 1 



(71) 



A;(r) = ;^^^Ml-^-;^ l - 241og 



6 288 



(72) 



(73) 



where dots stand for higher orders e^^^^^. 
G{r) = 27r^2-Fi ^) is a hypergeometric 



function and Li2(u) = ^ is a polylogarithmic 

function. R = (Qc/4)^/* is the overall radius of the 
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internal space (and of the AdS factor in the unflavored 
case). Notice that J enters always multiplied by in 
the equations of motion (51 )-(56 1, hence only the leading 



contribution in J is relevant in the solution. 



The horizon radius satisfies Gtt{f) 



(r — rh)bo{r) + 0{e\, S^) with 6o(r/i) ^ 0. As mentioned 
before, the radial coordinate r is derived after the 
solution has been obtained in cr, by demanding that 
Gxx{t) — / B? just as in the unflavored supersymmet- 
ric solution. Its relation to the coordinate a in which the 



equations of motion acquire the form given in (51|-(56) 
is also perturbative 



(1 



,1/4 



16 



1 + —euS^ { eK'^ - 2 - 2 tanh"^ 



4r^a(l + coth(2r^a))(l 



log 2) + log 4 



4rtcr 



(74) 



The one above must be understood as an effective 
IR solution, where the UV has been decoupled (i.e. all 
terms of the form r^/rg have been dropped). It should 
not be used for computing physical quantities at energies 
much higher than the plasma temperature. 

Let us conclude this section by summarizing the regime 
of validity of the solution. Apart from the usual bounds 
A^c ^ 1, A 3> 1 needed for the gravity approximation 
to be reliable, we must have 1 ^ N j <^ -\/A, where the 
lower bound comes from the requirement of validity of 
the smeared D7-branes configuration, while the upper 
bou nd is due to the use of the Abelian DBI (see sec- 
tion II E I . Finally, the perturbative solution is reliable if 
eii ^ 1, ^ ^ 1; the first condition allows to push the 
Landau pole far in the UV and have a predictive small 
temperature solution. Note that, as we will show in the 
next section, S ~ fi/VXT, so the requirement (5^1 
means that the chemical potential must not be para- 
metrically larger than the temperature T. 



V. THERMODYNAMICS 

In the previous section we have constructed a family 
of solutions consisting of black D3-branes and D7-branes 
smeared along the directions transverse to their world- 
volumes. These branes are dressed by a set of fields 
(scalar, vectorial and higher rank tensor forms), and 
depend parametrically on the coefficients eh and S. 
The size of the event horizon is dictated by these two 



parameters, and therefore thermodynamic properties 
associated to the horizon, for example the temperature 
or the entropy density, are given in terms of them. 



Given the solution we presented in (66)-(73), pertur- 



bative in Ch and S, we will perform the thermodynamic 
analysis up to order and S^, which suffices to catch 
effects due to the backreaction of the fundamental mat- 
ter. Consequently, all the results given in this section 
will have corrections of order 0(e^, <^^) which we will not 
write explicitly. Moreover, we will focus on the near-IR 
part of the geometry, equivalent to pushing the Landau 
pole to infinity (i.e. considering the limit — >■ oo as 
commented in the previous section). Thus, the results in 
this section have to be understood as providing correct 

answers up to subleading corrections O ^ 1. 

The position of the event horizon, r^, is determined 
by the solution to the equation b{rh) — 0. For finite 6 
there are two roots of the blackening function, and in 
this work we focus on the larger one. As the value of 6 is 
increased the two roots approach and, eventually, they 
are expected to merge, giving rise to an extremal horizon 
with zero temperature and finite entropy density. This 
situation is beyond the regime of validity of our solution, 
though, and will not be pursued. 

After Wick rotating time and imposing the absence 
of conical singularities in the geometry, we obtain 
the temperature (the inverse of the periodicity of the 
Euclidean time) 



■kR? 



1 - 



1 + 5^(2 -log 2) - 
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384^ 



-, 2(37 log 2 -26) 
13 



(75) 



The entropy, 5, is given by the Bekenstein-Hawking law and is proportional to the area of the horizon (an 
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eight-dimensional, fixed Euclidean time, fixed r = hy- 
persurface). As the solution considered here is an ex- 
tended black brane, this area is infinite due to the in- 
tegration along the spatial coordinates. To avoid the 
volume divergence we will consider intrinsic densities, 
s = S/ J dPx, and the (finite) entropy density reads 



2Vol(X5 



1 
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l + -e,{l + 5') + -4{l + 5') 



(76) 



where we have inverted ( 75 1 to express the position of the 
horizon in terms of the temperature and Vol(X5) is the 
volume of the five-dimensional Sasaki-Einstein manifold. 
The leading ©(e^J") term in expression (76) is the 



known unfiavored result. This will be always the case in 
the results that follow. The 0{e)^ 5^) term coincides with 
the result obtained in the quenched approximation.'^^ 
This agreement may seem surprising due to the different 
strategies followed to derive the result. In the present 
work the black hole size increases (with respect to the 
unfiavored setup) due to the presence of fundamental 
matter in the system. Furthermore, the branes de- 
scribing the fundamental matter are distributed along 
the directions transverse to the worldvolume of each 
individual D7-brane, explicitly breaking the symmetry 
group U{Nf) — ^ U{l)^f. On the contrary, in the 
quenched calculation a set of Nf <C Nc coincident probe 
D7-branes on the background given by the Nc D3-branes 
was considered .1221 Considering additivity of the system, 
the entropy associated to these degrees of freedom was 
calculated from the free energy and then added up to 
the entropy associated to the adjoint matter. A similar 
discussion holds for the 0{efi (5^) term and the results of 
the quenched approximation in the presence of a finite 
baryon density or a chemical potential.^'' Finally, terms 
that are higher order in or 5 are new results. 

The next quantity that can be derived just with the 
metric obtained in the previous section is the ADM en- 
ergy of the black hole. This quantity is not defined at 
the horizon but at the UV cutoff Ts — > oo. Formally it is 
given by 



Eadm — 



1 



dSabK 



ah 



(77) 



where Kab is the extrinsic curvature of a hypersurface 
at constant time and radius. This expression diverges 
in the — >■ oo limit and must be regularized. A use- 
ful subtraction scheme is given by considering Eadm — ^ 
Eadm — -E-adm where E^^Jyj^^j corresponds to the ADM 
energy of the solution with T = i5 = 0. Details of 
this calculation can be found elsewher^^ and we quote 
the final result for the subtracted energy density s = 
[EADM-Ef^,,)/ ^d^x: 



8Vol(X5) 



1 + 26' 



(78) 



Considering only 0{6^) terms in (78) the thermody- 
namic relation de = Tds is satisfied provided 



dT 



T 



(79) 



This expression can be derived independently from the 
fact that the gauge coupling constant = ingsc'^'' Nc 
is now scale dependent, as the dilaton <i>/i — $(r/i) runs 
at first order in e/i, see (71). This is dual to the log- 



arithmic running of the gauge coupling induced by the 
flavors, dXh/dT = e/jA^/T + ■ ■ ■ . Actually, the solution 
for the dilaton at this first order in is not dependent 



on 6, and therefore the relation (79) is valid when 0{S ) 



terms are considered. The inclusion of the 6 parameter in 
the discussion leads to the inclusion of a chemical poten- 
tial/charge density term in the thermodynamic relation, 
de = Tds + ^ duq . The determination of Uq , the quark 
density of the system, can be carried out holographically 



from the action ( 36 ) as 

_ SS _ 7r5/2Vol(A:3)7V/7VcV 
~ 5F^ ~ 16Vol(X5)3/2 



8 



(80) 

where we have used equation ( [75| ) and = Qc/'i- At 
leading order this expression coincides with the result 
obtained in the quenched approximation.t22l 

The product of the quark density, Uq, with the chemical 
potential, /i, can be obtained by considering the thermo- 
dynamic potentials in the canonical and grand-canonical 
ensembles, / = e — Ts and uj = e — Ts — fiUq respectively. 
These are related to each other by a Legendre transform, 
and with the Lagrangian (37) by an on-shell evaluation, 
as we will see later. Accordingly, one is able to obtain an 
expression for fj, Uq by evaluating on-shell the term 



/in. 



1 



5C 

Ja 



'- A' 



(81) 



and with ( 80 ) we can disentangle the expression for the 

r'(i+!) . 



chemical potential 

4V0l(X5)l/2 

(82) 

and once again the leading term was known from the 
calculation in the 't Hooft limit.'^Sl 

With this, the satisfaction of the relation de = Tds + 
jiduq at order 0{e\5^) is guaranteed provided 




6 
T 



3 



2 



(83) 



that we work at fixed Uq, such that the variation of 



(5 in ( 80 1 is compensated by the variation of the temper- 
ature and the gauge coupling. This in turn is a conse- 
quence of the fact that we are checking the thermody- 
namic relation in the canonical ensemble. To check it in 
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the grand-canonical one we would need to compare the 
ADM energy not to the T = 6 = case, but to the ex- 
tremal one in which b(rfi) — b' {rfi) = at fixed chemical 
potential /io- In this case the thermodynamic relation 
to satisfy is d{e — Sext) = Tds + — fio)dnq, however, 
this subtraction scheme is beyond the range of validity 
of our solution and we do not pursue it. As a result, 
in the grand-canonical ensemble the variation with the 
temperature at fixed chemical potential is obtained from 



equation ( 82 ) , resulting in 




S 
T 



eft, 
2 



-Oiel) 



(84) 



In the previous analysis we have used the definition 
of the free energy in the canonical and grand-canonical 
ensembles. The right-hand side of these definitions have 
been discussed in some extent now, and a relation of the 
left-hand side with the on-shell evaluation of the action 
was outlined. Next, we proceed to check that these two 
calculations agree with each other. 



The Lagrangian Lid m (39 1, depending on F123 ^ S, 



is the appropriate one for computing the free energy in 
the canonical ensemble, where the natural variable to 
consider is the quark density. The action derived from 
this Lagrangian has to be supplemented with a Gibbons- 
Hawking term to deal with a well-posed variational prob- 
lem. Despite the addition of this term, the evaluation of 
the on-shell euclidean action presents divergences in the 
UV {vs 00). These divergences can be tamed with the 
inclusion of an appropriate series of counterterms or by 
the subtraction of a suitable reference background. As 
we did in ( 78 ) for the ADM energy of the black hole, we 



subtract the setup with T = 6 = 0, such that 



F 



Jd^x 



— It — Io 



(85) 



where It = 'S'//s + S'//-|-5G//|on-shcii is the on-shell action 
(with the corresponding Gibbons-Hawking term) for the 
black hole solution and Iq is the corresponding quantity 
in the case with no black hole but periodic Euclidean 
time. The periodicity of the time in the last case can be 
arbitrary, however, we must ensure that the geometries 
we are comparing are matched at the UV cutoff (which we 
eventually send to infinity), and therefore the periodicity 
of the Euclidean time in the S — setup, f3o, is related 
to the periodicity of the black hole Euclidean time, /?, 
as (3o — \/b(rs)P. With this scheme the resulting free 
energy is finite^ and given by 



f^e-Ts 



(86) 



Using the expressions (79) and (84), —df/dT = s holds, 
resulting in a non-trivial cross-check. 

To study the grand-canonical ensemble we must con- 
sider the Legendre transform (with respect to A'^) version 



of (37). This operation is accounted for by just a change 
of sign in the -^123 term, and a similar discussion to the 
one performed before for the canonical ensemble leads to 



Ts 



II n 



1 ' 



(87) 



and —duj/dT = s. In this case, the variation of S with 



the temperature is determined by expression ( 83 ) 



Let us summarize what we have attained so far. Given 



the perturbative solution presented in (66)-(73l for the 



intersection of D3 and D7-branes in the Veneziano limit, 
we have succeeded in the calculation of the thermody- 
namic potentials in the canonical and grand-canonical 
ensembles by the use of the holographic dictionary. The 
explicit confirmation of the thermodynamic relations 
depends crucially on the variation of i5 with the tem- 
perature, which is an ensemble dependent property, 
and on the running of the gauge coupling due to the 
non-vanishing beta function that can be deduced from 



Once the potentials in the (grand-)canonical ensembles 
are known we can proceed to study the thermodynamic 
stability of the solution. A first check consists on deter- 
mining the determinant of the matrix of susceptibilities 



'^,T,T ^,iJ.,T 



where the comma stands for a partial derivative. This 
matrix is symmetric and — = x is the standard 
quark susceptibility. The components of the matrix are 



3^5 



UJ^T,T 



2Vol(X5) 

TrVoKATa) 
+ 4Vol(Xj 
7rVol(A:3) 



NfN, I? 



11 

24^ 



2Vol(X5) 
7rVol(A:3) 
4Vol(X5) 



NfN, fi T 



1 + jCh 

6 



6 



, (89) 
, (90) 
(91) 



and the determinant is parametrically positive, signaling 
thermodynamic stability. 

In appendix [A] it is shown that for A/" = 4 SYM 
(Vol(A3) = 2^,Vol(X5) = TT^) at zero coupling 
x(A = 0) = NfNcT'^, so that up to first order in we 
have ~^ co)/x(A = 0) = ^[1 + ge/i]- This means that 
the effect of the dynamical flavors is to increase the ratio 
of the susceptibilities of the interacting over free theory, 
pushing it above the value 1/2. The analogous situation 
in QCD is such that in the quenched case this ratio 
is close to one at relatively small temperatures above 

(see also the comments in section 3.1.1 of Ref. 
while in the unquenched case the ratio is reduced to 
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approximately 3/4 at T ^ 2Tc^ Thus, on one side the 
effects of dynamical flavors is qualitatively opposite in 
the two cases (flavored TV = 4 SYM and QCD); on the 
other side, they make the numerical values of the ratios 
in the two theories closer to each other, much like the 
corresponding ratios concerning other thermodynamic 
variables (s, e, uj). 

We can study also the change of the internal energy as 
we vary the temperature, described by the heat capacity 
of the system, cy = de/dT. The specific heat can be 
determined both at fixed chemical potential /i or quark 
density Uq and, for a thermodynamically stable phase, 
the result should be positive. Indeed, a direct evaluation 
considering ( 83 ) or ( 84 ) , depending on whether we fix rig 



or fj., gives the positive result 



2Vol(X5) 
2Vol(X5) 



1 



'All. 



h 

24 

2 

24 
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With the heat capacity at constant chemical potential 
we can read off the speed of sound 



1 



6 



(92) 



This result is parametrically positive and always lower 
than the conformal result = 1/3, in agreement with 
proposals in Ref. |42l Consequently, we see that confor- 
mality is broken at order by the addition of fundamen- 
tal matter. The fact that this effect is not seen at linear 
order in e/i is due to considering massless flavor degrees of 
freedom. Furthermore, the deviation from conformality 
is not affected by 5. To conflrni this point we can obtain 
the interaction energy of the system e — Zp (where p is 
the pressure) , which is a quantity that vanishes for a con- 
formal theory in three spatial directions. Using p = —lj 
we obtain 



e~ip _ n'^N^ 
T4 ~ 16Vol(X5 



(93) 



which yet again is proportional to and independent of 
6. 



the small parameter eh- In the following, focusing on 
the simpler uncharged case, we will review how in such a 
scenario all the hydrodynamic transport coefficients, up 
to second order in the hydrody namic al derivative expan- 
sion, can be explicitly evaluated.'^^IISlThe coefficients are 
given in terms of a single parameter (e.g. the speed of 
sound) and the results extend to any holographic plasma 
where conformal invariance is slightly broken by a small 
marginally (ir)relevant deformation dual to a single scalar 
field in the bulk. To these systems, which we will call 
"holographic marginal plasmas" we devote the following 
section. 



A. Holographic Marginal Plasmas 

Let us consider a 4d strongly coupled plasma dual to 
a simple 5d gravity action including just the metric and 
a scalar field if with some potential V{ip) 



(94) 



It can happen that the potential has an AdS^ (black 
hole) minimum at iy9 = const, (say, at ip — 0), i.e. that 
V{f) « -12 + (mV2)(/j2 _^ 0((^3) (setting the AdS 
radius to one). In this case we know that the scalar 
field is holographically dual to a 4d operator whose 
dimension A, around the conformal fixed point dual to 
the AdS background, is given by A(A — 4) = m^. If 
A < 4 (resp. A > 4), the operator is relevant (resp. 
irrelevant). If A = 4 the operator is exactly marginal 
and docs not drive the theory away from the conformal 
point. Examples of exactly marginal deformations are 
the "/3-deformations" .'^^ 

Marginally (ir)relevant operators, instead, break con- 
formal invariance, driving logarithmic flows of the cou- 
plings, despite their classical dimension being equal to 
four. The way to holographically account for such oper- 
ators in the simplifled model considered here is the fol- 
lowing: let us assume that the potential, for small values 
of some parameter 7, has a leading linear term 



V{ip) w -12(1 + -fip) + C'(7V^) 



(95) 



VI. HYDRODYNAMICS 

In the D3-D7 field theories with massless flavors which 
are the subject of the present work, conformal invariance 
is broken at the quantum level. The breaking is due to 
marginally irrelevant terms in the action, accounting for 
the coupling of the fundamental flavor flelds (introduced 
by means of the D7-branes) with adjoint (or bifundamen- 
tal) fields of the quiver theory on the D3-branes. In the 
previous sections we have shown how to treat the con- 
formality breaking effects in a perturbative expansion in 



At leading order the scalar field is thus massless and 
corresponds to an operator of classical dimension four. 
Moreover the 5d action does not admit an AdS mini- 
mum with (/p = 0, if 7 7^ 0. If I7I <C 1 we can easily 
find a perturbative solution for metric and scalar. At ze- 
roth order in 7 the metric is AdS (black hole) and, say, 
ip = 0. At first order, the equation of motion for the 
scalar field is solved (assuming ip to depend only on the 
radial variable of the background) by 

T 

(p{r) — —3-f log — . (96) 
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Here we have considered the non-extremal background 
with horizon radius r^- 

The logarithmic running of the scalar field and the 
AdS /CFT radius/energy relation r = E (which is 
still valid at leading order in 7) essentially realize the 
logarithmic RG flows expected in the 4d theory. If, as it 
happens in the simplest models where is the dilaton, 
the 4d running coupling is holographically given by 
gl^ = e"^, then models with 7 > (resp. 7 < 0) are in 
the marginally relevant (resp. irrelevant) class. 

There are many known string embeddings of the sim- 
ple model considered above, extending from examples 
with "fractional" brane^^ in type OB string theory'^ 
and in the type IIB conifold modeP^ - which sit in the 
marginally relevant class -, to the examples with flavor 
D7-branes considered in the present review, which are in 
the marginally irrelevant class. In these cases, with just 
one parameter driving the deformation (roughly speak- 
ing, the number of fractional or flavor D-branes added to 
an otherwise conformal model), the 5d gravity action ef- 
fectively reduces to the single scalar one considered above 
when the parameter is taken to be very small. This scalar 
is a string modulus dual to the marginally (ir) relevant 
term in the 4d field theory action. In general, of course, 
the 5d reduction of a type IIB gravity action gives rise 
to various scalar fields, but at leading order it is just 
this scalar which plays an active role. The other ones, 
at leading order, are just "frozen" at their zeroth order 
values. For the D3-D7 uncharged plasmaP^jl and the cas- 
cading conifold one,'^ this has been explicitly shown in 
Ref. [12] We report the analysis for the uncharged D3- 
D7 plasma in appendix|B] from which it follows that, at 
leading order. 



so that the simple equation of state 



(/3 = $ - , 
where <& is the dilaton. 



eh 

' 3 ' 



(97) 



B. Equilibrium Properties: Equation of State 

One of the simplest way to realize a scalar potential 
going linearly at leading order is by choosing it to have 
an exponential form, V{ip) = — 126^*^, as it happens 
in the Chamblin-Reall models.!^ Of course, the details 
of the non-leading terms in 7 are irrelevant if one 
wants to focus on a leading order analysis. Anyway, 
the Chamblin-Reall models have a first advantage of 
being exactly solvable and so, in particular, an exact 
black-hole solution is known to all orders in 7. The 
thermodynamics of the dual 4d models thus easily 
follows. In particular, as it has been discussed in Ref. 
H51 the speed of sound is given by 



P = 



T 
2 



(99) 



dp 

de 
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follows; here p is the pressure and e is the energy density. 
These expressions explicitly show how the parameter 
7 weighs the conformality breaking effects: for 7 = 
we get back the conformal result = 1/3, i.e. e — 3p. 
Indeed, using 7 = — e;i/3, we see that the speed of sound 
precisely matches that obtained for D3-D7 plasmas with 
Eh <C 1 as discussed in the previous sections. 

As has been discussed in Refs. i49i and iSO^ Chamblin- 
Reall models in d + 1 dimensions have a further nice 
property: for particular values of the coefficient of the 
exponential in the potential, they can be obtained from 
dimensional reduction on a 2(t — d torus of Einstein 
actions with negative cosmological constant in 2(t -I- 1 
dimensions. This happens when the parameter ct, which 
determines the coefficient in the exponential together 
with d, is half-integer. For these values of <t, one can then 
start from the well-known AdS^a+i solution and its dual 
energy-momentum tensor (at or near equilibrium) and 
obtain the energy-momentum tensor for the dual to the 
Chamblin-Reall model by simple toroidal dimensional 
reduction. 

For instance, let us just focus on the equilibrium 
properties in the d = 4 case. The higher dimensional 
action has an AdS2a+i (black hole) minimum, dual to 
a 2(j-dimensional (thermal) CFT. The equation of state 
for the latter is just given by e = {2a — l)p as it follows 
from the vanishing trace of the stress-energy tensor. 
The toroidal dimensional reduction does not affect the 
equation of state which will thus be the same for the 
4d plasma dual to the Chamblin-Reall model. Using 
the previous result, we see that 2a — 1 = 6/(2 — 37^). 
Thus, if 7 ^ 1, cr « 2 -I- (9/4)7^. As we have previously 
recalled, the analysis requires cr to be a semi-integer and 
this is certainly not realized for any 7. 

However, a crucial observation in Ref. [SUI is that, 
from the point of view of the theory in d -I- 1 dimensions, 
the equations are smooth in the parameter a > d/2 
(at a — d/2 the action is singulaJ^. Then one can 
proceed as follows. One starts from a Chamblin-Reall 
model in d + 1 dimensions for whatever a > d/2 and 
performs the continuation (which is smooth) to the 
nearest value a which is semi-integer. The latter theory 
is the compactification of a theory admitting a AdS2d+i 
solution, so its dual energy-momentum tensor, which will 
be a function of a, can be calculated straightforwardly. 
This energy-momentum tensor can thus be continued 
(smoothly) back to the one of the theory corresponding 
to the original value a. For the case d = 4 and 7 <C 1 we 
can thus choose a = 5/2. 
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Notice that this procedure ahows to get the full stress 
energy tensor at or near equilibrium. In the latter case, 
we can thus deduce the full transport coefficients of the 
4d plasma dual to the 5d Chamblin-Reall model, just by 
dimensional reduction of the hydrodynamic stress tensor 
of an holographic higher dimensional conformal plasma. 



SYM. All the others coefficients, i.e. the bulk viscosity 
C and the second order coefficients K*,r*, A4,Tn ("bulk" 
relaxation time), i^i, Csj '?4; Cei are only defined in 
non-conformal plasmas. 

D. Hydrodynamics of holographic marginal plasmas 



C. Hydrodynamics: Generalities 

Hydrodynamics is the effective theory of long wave- 
length, low frequency fluctuations around local thermo- 
dynamic equilibrium. In the uncharged case relativistic 
hydrodynamics is determined, up to second order in the 
derivative expansion, by seventeen transport coefficients, 
fifteen of which are possibly independent .l^^^t^ On a gen- 
eral space with metric gi^j^, the energy momentum tensor 



(100) 



where A'"' = g^'^ + u^u^ , is determined by the energy 
density e, fluid velocity u'' {u^u^ = —1), the transport 
coefficients in its "viscous shear" part: 



V • u 



. V • M 



+t]t; — -ff'"' + A4V<^ log sV''> log s , 
and in its "viscous bulk" part: 



(101) 



n = -C(v • u) + CmDiw ■ u) + a^^'fTp^ + 6(V • uf 

Hs^^'^l^. + U'^i log sV^ log s + 



<5R + i&u'^uf^R. 



'a/3 ■ 



(102) 



while the pressure is given by the equation of state p{e). 
The various structures in ( 101 1 and ( 102 ), apart from the 



obvious Riemann and Ricci tensors and scalar curvature 
{R^"'P'^,R^"',R), are given by 



D = u^V. , 



V^u^- -A^'^(V-u), 
3 



(103) 



and for a generic tensor A^" the notation was used that 

{A^"') = A<>"'> = lA''°'A''^{Aa,p+Apa)-\A'"'A''^Aap . 

(104) 

Finally, s is the entropy density, while — dp/de is the 
square of the speed of sound. 

The shear viscosity rj and the second order coefficients 
("shear" relaxation time), k, Ai,A2,A3 are the only 
ones defined in conformal fluids, such as the one of A/" = 4 



Just like any theory with a dual gravity description 
satisfying the general requirements established in Ref. 
[S31 4d marginal plasmas dual to (small 7) 5d Chamblin- 
Reall models have a shear viscosity over entropy density 
ratio given by (in units h — Kb — 1) 



1 

47r 



Moreover, according to the results in Ref. 
els have a bulk viscosity given by 
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(105) 
these mod- 

(106) 



so that in the D3-D7 case, ^ = at leading order (in 
appendix|B]we will cross-check this result using a general 
holographic formula for the bulk viscosity proposed in 
Ref. I55() . Comparing this expression with that of the 
speed of sound in cq. ( 98 1 we see that the bulk viscosity 
saturates the bound 



c 



> 2 



1 



d-1 



(107) 



proposed in Ref. 132! for any d + 1 dimensional plasma 
with gravity dual. As observed in Refs. fT^ and fT3, this 
is a generic feature of holographic marginal plasmas at 
leading order. 

In order to compute the second order hydrodynamic 
transport coefficients, one can take advantage from 
the fact that, as previously discussed, for 7 ^ 1, 
5d Chamblin-Reall solutions can be obtained from 
reduction and analytic continuation (from a = 5/2 to 
a ~ 2 + (9/4)7^) of an AdS2a+i background, whose dual 
conformal hydrodynamics was considered in Ref. '56' Let 
us skip the details of the (simple) derivation, for which 
we refer to the original papeif^^ (see also Refs. [53land [57|) . 

The final result is shown in table |T] where we give 
the full list of transport coefficients for a holographic 
marginal plasma, at leading order in the conformality 
breaking parameter 



(108) 



In particular to get the transport coefficients of the un- 
charged flavored D3-D7 plasmas, to second order in e/j, 
it suffices to notice that 



6 



(109) 
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According to what has been previously commented (see to the cascading plasma as wellj^ at leading order in 
also Ref. 12 j, the results in table |l] automatically apply Set- 
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TABLE I. The transport coefEcients, in the notation of (100l-(102l, for a marginally (ir)relevant deformation of a conformal 
theory, at leading order in the deformation parameter Set = (1 — 3cJ). The holographic equation of state is £ = 3(1 + Scb)p- For 
D3-D7 plasmas Set 



4/6. 



A comment has to be devoted to the two relaxation 
times Tt^jTyi- At leading order in the conformality 
breaking, the bulk relaxation time Tn is not proportional 
to the bulk viscosity. The shear relaxation time r^r, 
moreover, depends on the speed of sound. Interestingly, 
the magnitude of the shear relaxation time is larger 
than the one of the bulk relaxation time rn; moreover, 
if we naively extrapolate the qualitative temperature 
dependence of the QCD speed of sound to the present 
formulas, r^r is found to increase more steeply than rn 
as the temperature is reduced, contrary to the common 
speculations about the behavior of these two coefficients. 
The discrepancy could depend on the strong coupling 
regime we are considering. 

In addition, using the above results it is easy to verify 
that the relation 



4Al + A2 = 2l]Ty, , 



(110) 



holds, at first order in Srh- It has been shown in Refs. 



and i59^ that (110) is satisfied in all the known examples 
of conformal plasmas (in d > A spacetime dimensions, 
with or without conserved global charges) with a dual 
gravity description. Our results provide a unique validity 



check of the above relation in non-conformal settings. 

As an independent check of the results shown in table 
[ij we will provide, in appendix [C] a direct holographic 
computation of some (combinations of) the first and 
second order transport coefficients. 

Let us conclude by mentioning that in Ref. [131 holo- 
graphic marginal plasmas have been used as toy models 
for the initial stages of the evolution of the Quark-Gluon 
Plasma produced in heavy ion collisions at RHIC and 
LHC. Indeed, if the temperature of the plasma is some- 
what higher than the temperature for the QCD crossover 
between confined and deconfined phases, it is sensible to 
approximate the QGP as a slightly deformed conformal 
fluid. The results in table|l]thus give all the transport co- 
efficients in the toy model, once only the speed of sound 
at some temperature is given (from lattice data). 



E. Negative Refraction and Additional Liglit Waves 

Thermodynamic quantities and hydrodynamic trans- 
port coefficients are sufficient to determine quantitatively 
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some interesting optical properties of the D3-D7 plas- 
mas. With "optical properties" we mean the response 
of the system to an external U{l)em electromagnetic 
incident wave. At leading order in the electromagnetic 
coupling g, the C/(l)em can be treated as global and 
we can use our charged solution as a model of an 
electrically charged, strongly coupled relativistic plasma, 
i.e. we treat U{1)b as the U{l)em-, and the fundamental 
matter as ?7(l)em-charged particles. With respect to the 
R-symmetry charged solutions, the model at hand has 
the advantage that, as in ordinary "phenomenological" 
systems, only a fraction of the fluid constituents is 
charged. 

As a first optical property, it has been shown in Refs. 
1501 and EH that any charged plasma exhibits negative re- 
fraction of light for small enough frequencies (holographic 
studies also appear in Refs. fHllS^MS)) . This means that 
the phase velocity and the energy flux have opposite di- 
rections, resulting in a number of non-standard optical 
properties. The critical frequency for the negative re- 
fraction is (we consider here only the A/" = 4 SYM case) 



Re[n] 



P 



1 



(111) 



Since riq/T^ ^ 0{NfNc^/XfiS), the window of frequen- 
cies where negative refraction is present is 0{NjXhd'^). 

Moreover, it is known that the same phenomenon that 
gives rise to negative refraction, i.e. spatial dispersion, 
is responsible for the possible presence of so-called "Ad- 
ditional Light Waves" (ALW) in certain materials: for 
each incident light wave, there could be more than one 
propagating waves, differing for their refractive index. 
This phenomenon has been studied holographically in 
Ref. 66 in the R-symmetry charged solutions. Again, 
in the plasma regime the phenomenon is (obviously) 
entirely determined by thermodynamics and transport 
coeflticients. 

For the solution at hand, the resulting effect is very 
similar to the R-symmetry charged case: using formulas 
in Ref. |66l one can check the presence of at least one 
ALW besides the normal wave. In figure |4] we report the 
plot of the real part of the refractive index Re[n] as a 
function of the frequency w for the two waves. The bot- 
tom line in the plot refers to the ordinary wave, which 
exhibits negative refraction (signaled here by a negative 
Re[n]) for frequencies smaller than the critical one given 
in ( 111[ ). The upper line refers to the ALW, with positive 
refraction. The middle line is an "effective refractive in- 
dex" measuring the relative magnitude of the two waves: 
its interpolating behavior between the two other lines sig- 
nals the fact that in the intermediate frequency regime 
both waves are propagating with comparable amplitudes. 
Unfortunately, as in all the holographic models analyzed 
so far, the imaginary part of the refractive index is very 




FIG. 4. Real part of the refractive index, as a function of the 
frequency, for the ALW (top line) and the normal wave (bot- 
tom line); the middle line is the "effective refractive index". 
In this plot we used the parameters Nc = 100, Nf = 20, = 
50,g = 1,(5 = 0.1. 



large in the interesting frequency regime, producing a big 
dumping of the propagating waves. 
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Appendix A: The U{1)b Susceptibility of the Free Flavored 
A/" = 4 SYM Plasma 



Let us consider the M = 2 theory obtained from 
TV = 4 SU{Nc) SYM coupled with Nf massless funda- 
mental hypermultiplets. The fields in the latter have the 
same charge q = 1 under the baryonic U{1)b- All the 
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other adjoint fields in the theory are naturally uncharged. 
In order to compute the total U{1)b susceptibilitjE3 

duB 



X 



(at^ = 0) , 



(Al) 



for the free theory, we just need the expression for the 
net total baryon charge density ub- For each field the 
net charge density is given as the difference between the 
particle and the antiparticle density. Thus, for complex 
bosons of charge q under a [/(I) 



Tib 



(Pk 
(27r)3 



[Fk{Ek ~ qii) - FtiEk + q^i)] , (A2) 



and for Weyl fermions (note: for Dirac fermions there is 
a further overall factor of two) 



(27r)3 



[Ff{Ek-q^i)-Ff{Ek+q^l)] . (A3) 



In the expressions above E^ = \k\ = k since we want to 
consider massless particles. Moreover 



Fbix) 



1 



1 



3/3=^ + 1 ' 



(A4) 



are the standard Bose and Fermi distributions with 
/? - l/T. 

From the expressions above it is easy to compute the 



susceptibilities using the definition (Al). The results are 



Xb^q 



Xf 



•■II 
6 



(A5) 



In our case, counting the fields charged under U{1)b, we 
have 2NcNf complex scalars and 2NcNf Weyl fermions. 
Thus the free baryon number susceptibility is ((/ = 1 
here) 



1 1 



Xi\ = 0) = 2iV,iV/ + = N.NfT 



(A6) 



In the planar infinite 't Hooft coupling limit the same 
theory, in the quenched approximation, has (see Ref. 



X{X ^ oo) = -N.NjT' 



(A7) 



so that x{X ~^ co)/x(A = 0) = 1/2 + . . . as it happens 
for the R-charged case (see for example appendix A in 
Ref. _69 and the comments in Ref. ■ 

Just for comparison, let us consider the free SU{Nc) 
QCD with Nf fundamental Dirac fermions. The to- 
tal U{1)b susceptibility of the related Stefan-Boltzmann 
quark-gluon gas is readily evaluated as 



XQCDiSB)^N,Nf2 



1 



':NrNfT 



(A8) 



Appendix B: The 5d Effective Action in the Uncharged 
Case 



The ten dimensional action (20), (21) can be reduced 



to a five dimensional effective action with a suitable 
ansatz for the various fields. We will only discuss the 
uncharged case, so that F = i^(3) = B = Q. Moreover, 
-^"(1) = Qf{dT + Ake) and i^(5) = Qc{l + *)e5 (the last 
symbol is of course the volume form of the internal mani- 
fold) . The reduction ansatz of the metric is of the forrrP^l 



2 

KE 



+e 



2(4™-/) 



(dr + Ake) 



(Bl) 



We will use Greek indices v for the five dimensional 
space. In the parametrization above 



/ = -ilog(/iJ5*i^ 



5 



(B2) 



and 



dsl 



(hiS^F)i 



-h-^'%de 



h-^l^dxl 



-h^'^S^'F^da^ 



(B3) 



Plugging the above ansatze in the ten dimensional ac- 
tion and integrating in the five internal direction, we ob- 
tain the five dimensional effective actiorP^I 



^5 



V0l(X5 



^ / d^x^J- det g 



40 



2Q{dwf 
(B4) 



where 

y($, /,«;)= 4e^^'+2-(eio-- 6- 

We have set Racls = 1 and a' = 1 for simplicity. In 
these units Qc = 4. 



, 40 f 

-6 3-'. 



(B5) 



Let us study the above 5d action perturbatively in 
th, writing each of the scalar fields as well as the 5d 
metric components as 4* = 4*0 + e^^i + At order 
zero the action has an AdS (black hole) minimum at 
f = w = and $ = const. = ^h- The field / is dual 
to an irrelevant operator of dimension A = 8 and the 
field w is dual to a vev for an irrelevant operator of 
dimension A = 6. The dilaton $ is dual to the insertion 
of a marginally irrelevant operator, actually the fiavor 
term in the field theory (T = 0) superpotential.'^® This 
is the source term which is responsible for the breaking 
of conformal invariance at the quantum level. 

To better understand the role played by the various 
scalars, let us consider the quantity 



(B6) 
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in three different cases, where we identify with one of 



the three fields entering the potential (B5), taking the 
other two scalars fixed to their background values. At 
order we get 



contributes to the bulk viscosity, and using eq. ( B8 ) we 
find C/?7 
lYIl 



e^/9, confirming the result found in section 



V 



(B7) 



This indicated^ that, among the three scalar fields in 
the action, only the dilaton - i.e precisely the field dual 
to the source for the (marginally irrelevant) deformation 
driving our theories away from conformality - plays the 
role of "active" field in the game at leading order. The 
other two scalars, / and w, do not contribute at leading 
order and they can be fixed to their background values. 

All in all, at leading order in the conformality breaking 
parameter, the system is effectively captured by a simple 
5d dual gravity model with a single scalar field with 
linear potential. Actually, since is a constant up to 
higher order terms, the effective 5d action, at leading 
order, is mapped (see section VI ) into a Chamblin-Reall 



model^ with constant 7 = — ^ (the sign is chosen in 
order for the solution to the scalar field equation to 
match the one found for the dilaton in the main body of 
the paper). 

Using the holographic field/operator map found in Ref. 
[701 in Ref. |T2]it was shown that analogous considerations 
hold for the cascading plasma too, at leading order in the 
conformality breaking parameter. 



Appendix C: Alternative Computation of Transport 
Coefficients 



In this appendix we present the direct computation of 
some (combinations of) first and second order transport 
coefficient s, provid ing an independent check of the results 
in section |VI [ ^^I^'^ I Specifically, the pressure p, the trans- 
port coefficients i], k and a combination of the "shear" 
relaxation time r^r and k* are derived from the retard ed 
correlator of the tensorial fiuctuation of the fluicP^IIl 



= p-i'quj+ (t]t^ - I + k*) _ _^ O{q\oj^) . 

(CI) 

In this expression w is the frequency of the fiuctuation 
and q its momentum. The shear viscosity j] can be ex- 
tracted also from the dispersion relation of the vectorial 
(hydrodynamic) fluctuation mode 



CO = -I ^q^ + 0{q') . 



(C2) 



Moreover, the speed of sound c^, the bulk viscosity ( 
and a combination of the "shear" and "bulk" relaxation 
times , Tn are derived from the dispersion relation of 
the scalar (hydrodynamic) mode (sound channel) 



1. The Bulk Viscosity From a General Holographic 
Formula 

As a further check of the previous conclusions, which 
in section |VI| have been heavily used to study the hydro- 
dynamics of the D3-D7 plasmas, let us consider a general 
holographic formula for the bulk viscosity 



c 

V 



E 



ds 



dp"' 



(B8) 



found in Refs. |55] and [TT] In this expression the 
suffix h means "evaluated at the horizon", (j>i are the 
scalar fields with canonically normalized kinetic term 
in the 5d gravity action, and is the charge density 
holographically related to the a — th U(l) field in the 
5d action. In our case we have just the U{1)b charge. 
Finally, s is the entropy density. 

Let us just consider the uncharged case for simplic- 
ity. The scalar fields in the 5d action are ^,f,w. 
Their background values are $ ~ <i>/i -I- e/j log(r'/r/i), 
/ « — e/i/40, w « — e/j/60 at leading order. Now 
sd^/ds fa {rh/3){d^/drh) ~ e^/B up to subleading 
terms. The analogous quantity for / and w is of or- 
der . From this we readily find that only the dilaton 



io^Csq- iTq' + £ (c^r^// -l^^ + Qiq^), (C3) 
where 

The fluid modes are dual to fluctuations of the five 
dimensional gravity fields above the solution of section 
|IV| in the zero charge case. The action for the fluctua- 
tions is the one in formula (B4). In the calculations we 



systematically discard the contributions coming from the 
UV completion of the field theory, which are power-like 
terms suppressed as (r^/rg)". 



1. Fluctuations 

Following the standard procedure,'^ we assume that 
the perturbations take a planar wave form in Minkowski 
space, ^^(x^jr) = e-i(^t-<i^)\ij(^r). The perturbations 
can be classified according to their transformation under 
the little group 5*0(2), which is the remaining symmetry 
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of the system (rotations in the x — y plane) . We define 
,5<7„(x^r) = e-'("*-«-).9„(r)H«(r), 



r) = e-'("*-«-)5,,(r)H„„(r) , (m, n) ^ {t, t) 



-i{ujt—qz) 



-i{ujt—qz) 
-i{ujt—q 



B{r) 
'Mr)- 



(C5) 



We work in the gauge Hmiir) — 0. The system of per- 
turbations above includes a tensorial mode (Jixy), vecto- 
rial modes {Htx, T~Lzx, T^ty, T^zy) and scalar modes CHtt, 



n 



±± 



'^yy: 'Hzz, "Htzi C). Each kind of 



perturbation can be expressed in term of gauge invariant 
quantities under the residual gauge symmetry^'^ 

Tensorial 



Zt 
Zv 



qHtx 



Vectorial 

Scalar Zs ^ 2H 



^-Utz- 



1 - 



g'xx 



2 

, „ 9 fftt T, 

gxx 



B 



log [g, 
f 



■Ha 



Zc^C 



log' [g 
w' 



log' [5: 



■n 



(C6) 



The differential equations for all these fluctuations at the 
horizon admit solutions behaving as — , where 

ro = uj/{2rii) and = ttTq (Tq denotes the temperature 
of the unflavored A/" = 4 SYM theory). The index with 
negative sign corresponds to incoming wave boundary 
conditions at the black hole horizon, which in turn give 
the causal solutions we are interested in. We impose that 
the fluctuations vanish at the UV cutoff scale related to 
rg'. the results turn out to be independent of rg up to 
suppressed terms in powers of rh/rg. 



2. Tensorial Perturbation 

We scale n> Xhyd^, q A/^y^q, where q = q/{2rh) 
and Xhyd is a parameter keeping track of the order of the 
hydrodynamic expansion. Define 



Z. 



„4 

r4 



"^0 2 2 



j=0 k=0 



(C7) 



where higher order terms in and Xhyd, which we will 
not study, are not taken into account. From the action 



( B4 1 we derive the following equation for the perturba- 



tion 



ZrJ^ 



1 



Z^log'( 



gttgxx 



g-rr / 2 2 9tt \ „ „ 

I a; — q Zt — . 

gtt ^ gxx'' 



(C8) 



With the ansatz (C7| one can check that the only 



non-zero term in the solution normalized to one at the 
horizon up to first order in X^yd is 2!^^ = 1. At second 
order in X^yd the solution is too lengthy to be reported 
here but straightforward to obtain. 

Dirichlet boundary conditions cannot be imposed on 
the solution in the UV, showing the absence of a hy- 
drodynamic mode in this channel.l^ We can nevertheless 
write the hydrodynamic expansion of the retarded cor- 
relator once we evaluate the action on-shell. The latter 
is singular when evaluated at r = — >■ oo. To cure this 
divergence we have to add the following countertermJ^H 



Sbulk Sbulk + 



Vol(X5) 



d*£.2y/^K 



(C9) 



Vol(X5) 
2k\q 



-7(W[0]-ic[0]i?[7] 



where K is the scalar associated to the extrinsic curva- 
ture of the cut-off hypersurface, 7 is the four-dimensional 
metric at the boundary, C[0] is a function of the scalars 
and Wlcj)] is the superpotential 



=5/ 



Qf 



4to-|-<I> 



4e 



f+6w 



(CIO) 



from which the potential terms in the action ( B4 ) can be 
derived as 



80 V df 



1_ f dW 
40 V dw 



fdW\ 



1, 



(Cll) 

The function C[(/)] satisfies the differential equatiorP^l 

2^ l[ml)TW ^ ^~d^dw ^ ~d¥d^\~^^ ~ ' 

(C12) 

The exact form of C[0] is actually not needed in order 
to extract physical results, since the divergence balanced 
by this function goes in the UV as r^, being the next-to- 



leading order suppressed as 



it does not affect 



the finite part from which the hydrodynamic transport 
coefficients are obtained. The leading behavior, needed 
to cancel the divergence, reads 



23 

los' 



371 
23328* 



+ 0{r-^). (C13) 



The Fourier transformed, quadratic-in-fiuctuations, 
on-shell boundary action is 



5 = 



(CM) 



with Hk the boundary value of the fluctuation {J- is obvi- 
ously not to be confused with the D7-brane world-volume 
gauge field strength) . The retarded correlator of the cor- 
responding components of the energy momentum tensor 



G 



xy,xy _ 



I / dt(i'xe'^'^'-'''^Q{t){[Txy{t,x),Txy{QM) , 

(C15) 



24 



is related to the on-shell action bjl^ 

-2Im[J"(fc,rs 



^ixy,xy 



(C16) 



Plugging the solution of the equation of motion ( C8 ) 



for Zt, in the finite action (C9|, it is straightforward to 
derive the flux J^{k,rs), and so the correlator G^'^^ as 



Gr ''' = ^yf4l([l - 2.tr - 2q2 + 2w\l log2)] 



8Vol(X5; 
in> + 2q2 - 211)2(1 - log 2) 



-Eh 



(C17) 



24 + 19itt)-4q2 + 2tt)2(2 + 37r2 + 221og2) ^ 
192 



Comparing to (CI ) and using the expression for the tem- 
perature in ( |75[ ) (at zero charge), we confirm the value 
of the pressure p in section [v] (a t zero charge) and the 
results in Table [l] of section VI D| for the transport coef- 
ficients T], r^, K, K* . 



Vectorial Perturbation 



where 'yk,tk are the order coefficients of the dimen- 
sionless combinations 



7 = ttTo r , 



t = 



(C21) 



The relevant equations for the perturbations are 



f. <?' g'tt \ rr 



EOM 



f2^^i/, 
gxx 

1 77 ttEOM 

• EOM 



EOM 



^ gxx^ 

1 9tt " 



lo^axx^H^OM 

9tt 



0, 



(C22) 



log' glx 



H, 



EOM 



^9xx[g'xx'l>Bl'^ + 't>'B{^/9x 



1 9xx 9xx \J 9xx ) 



X H, 



EOM 
rt 



19tt jjEOM 

^gxx 



= 0. 



(C23) 



7" 



The equation in this channel reads 
'^1-log'f^ 

9tt 



log'(^)-log'(^)( 

V Qtt Qrr ' V Qtt I V 



1 



9u9. 

2 9tt 



^9u_^ 
^2 



[ijJ - q ]Zv = 0- 

9tt ^ gxx'' 
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For the vectorial fluctuation we can solve order by or- 
der with the scaling W — A^j^^tr, q — > Xhyd<^, imposing 
regularity at the horizon. The result is 



Zv = Cv 



1 -i- 



1 - 



„4 
'_h 
^4 



(C19) 



(1 + e, +e2) +0{rv,q^) 



From Dirichlet conditions at the boundary — oo we 
can read off the shear viscosity from the dispersion rela- 
tion (C2). This calculation is summarized in the mem- 



brane paradigm formula given in Ref. I77L and it gives 
the well-known ratio rj/s = l/(47r) with corrections in 
powers of r-^jr^ — 0. 



4. Scalar Perturbations 

We write for each scalar perturbation ZA=s.B,C,ip ths 
ansatz 



In these expressions, represents each of the scalars 
/, w,<I> (the form of their equation is the same) and (\)-b 
their background value. Moreover, the notation 0-^'-'*^ 
(and H^J^^^ and so on) stands for the corresponding 
equation for the scalar (and the fluctuation H^z and so 
on) in section 3 of Ref. [731 For the coefficient S we have 



Vu;2[q2(^4 „3^4) +3^4^2] 

4q^r*{q^ - oj^){r* - rl) 

Vw2[q2(r4_ 3^.4) ^3^,4^2] 



(C24) 



The solutions for the non-zero fluctuations entering the 
sound channel, satisfying the normalization at the hori- 
zon and Dirichlet conditions at the boundary in the case 
of Zg, read 



1 



^0,2 



\ogp 



12(1 -p4)' 
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7r2(p4-l) + 24(p4_i)log2p 



144^/3(1 - p4) 
-121ogp(4+(p4-l)log(l + zp) + 

ip' - 1) log [^{^ + p){p' 1)]) - 3{p' ~ l)Li2{p') 



^c,,feejq-2i X!7fce^q'+4X!^fc ^* 1^ 



k=0 k=0 
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fe=0 



Za = Ca\1 



4 \ ~ 2 2 



Y,T.^AU'et. (C20) 

j=0 fc=0 



where p = r/rh- We do not report the expressions for all 
the q^ coefficients of the solutions because of their very 
lengthy form. 

The Dirichlet conditions at the boundary give the dis- 
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persion relations in the first line of ( C20 ) with 



Cs,0 



70 



to 
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6 ' 



71 



3-21og2 



Cs,l = , 
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24^3 ' 96%/3 ' 
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which confirm the result for the speed of sound ( 92 1 found 



with the thermodynamics and, using ( C3 ) , ( 75 ) and ( 79 1 , 
confirm the results reported in Table |Ij of section VI D for 
the transport coefficients Ci ''"tt, '''n- 
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